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Abstract. We discuss high resolution VLT/UVES observations (FWHM ~ 6 km s" 1 ) from October 2002 (day ~ 5700 past 
explosion) of the shock interaction of SN 1987A and its circumstellar ring. A nebular analysis of the narrow lines from the 
unshocked gas indicates gas densities of (~ 1.5 - 5.0) x 10 3 cm" 3 and temperatures of ~ 6.5 x 10 3 - 2.4 x 10 4 K. This is 
consistent with the thermal widths of the lines. From the shocked component we observe a large range of ionization stages from 
neutral lines to [Fe XIV]. From a nebular analysis we find that the density in the low ionization region is 4 x 10 6 - 10 7 cm" 3 . 
There is a clear difference in the high velocity extension of the low ionization lines and that of lines from [Fe X-XIV], with 
the latter extending up to ~ -390 km s" 1 in the blue wing for [Fe XIV], while the low ionization lines extend to typically 

260 km s" 1 . For Fto a faint extension up to 450 km s" 1 can be seen probably arising from a small fraction of shocked 

high density clumps. We discuss these observations in the context of radiative shock models, which are qualitatively consistent 
with the observations. A fraction of the high ionization lines may originate in gas which has yet not had time to cool down, 
explaining the difference in width between the low and high ionization lines. The maximum shock velocities seen in the optical 
lines are ~ 510 km s" 1 . We expect the maximum width of especially the low ionization lines to increase with time. 
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1. Introduction 

The core-collapse supernova (SN) 1987 A in the Large 
Magellanic Cloud (LMC) is much closer to us than the SNe we 
normally observe. Detailed studies of SN 1987 A have therefore 
contributed tremendously to our understanding of this type of 
SNe. Its interaction with circumstellar gas has also c ontributed 
to our underst anding of Shockwave physics (see e.g. lChevalietl 
1997tlMcCravll2005h . 

The radioactive debris of SN 1987 A is at the center of a 
triple ring system where the two outerm ost rings are about 



three times the size of the innermost ring dWang & W ampler 



19921) . The rings are elliptical in shape with the inner ring cen- 
tered around the SN and the two outer rings are centered to 
the north and south of the SN forming a possible hour-glass 
structure. Since the inner ring (henceforth called the equato- 
rial ring, or jus t ER) appears to be i ntrinsically close to cir - 



cular in shape (IGould & Uza 1 1 1 998t ISugerman et al. 2005) 



Send offprint requests to: P. Groningsson 

* Based on observations collected at the European Southern 
Observatory, Chile (ESO programme 70.D-0379). 



its observed ellipticity is interpreted as a tilt angle to the 



line of sight of ~ 44° dPlait et al. II 19951: Burrows et al. 111995 



Suge rman et al. 112.0021) . How all the rings were formed is still 
unknown, but it seems clear that the ER is the result of the fi- 
nal blue supergiant wind interacting with previously emitted 
gas in the form of an asymmetric wind-like structure (e.g ., 
Blondin & Lundqvist 1 1 19931: IChevalier & Dwarkadasl |l995D . 



In the models by iMorris & Podsiadlowski I (120051) the strong 
asymmetry is explained in a merger scenario being due to a 
rotationally enforced outflow. In these models also the outer 
rings may be explained by the merger, and the rings are likely 
to have been formed during the last ~20,000 years before the 
explosion. Many details in both the observations and models 
are, however, uncertain. 

The observed ER was photoionized to a high state of 
ionization by the SN flash accompanying the SN breakout 



dFransson & Lundqvist Il989l) . t has since the breakout been 



cooling and recombin ing, giving rise to a mu l titude of narrow 



emission lines (e.g., iLundqvist & Fransson 



Il996l) . Already 
early it was realized that the ER should be reionized ~ 10-20 
years after the explosion when the expanding SN ejecta would 
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start to interact with the ER (Luo & McCray 1991; Luo et al. 



19941) . This event would cause a rebrightening of the ring and 
the first indication of this interaction became visible as a bright 
spot ( Spot 1) in 1997 on the north side of the ring at PA 
= 29° dGarnavich. Kirshner & Challis lll997tlSonneborn et al. 



1998). 



Modelling of Spot 1 (Mich ael et all [2000) suggests that 



it is an inward protrusion of the dense (~ 10 cm , 
Lundqvist & Fransson 19961) ER and that it is embedded in 



the lower density (~ 10 cm , Ch evalier & D warkadas 1 11995 



iity (~ l 

Lundqvist II 19991) H II region interior to the ring. The shocked 
gas drives a forward shock into this H II region and the bright 
Spot 1 is the result of the interaction of the shock and an in- 
ward protrusion of the inner circumstellar ring. When this hap- 
pened, slower radiative shocks were transmitted into the protru- 
sion and the radiation from these shocks appeared as a bright 
spot. 

The velocities of the transmitted shocks are given by V s ~ 
Vft(nHiiMsp) 3 , where Vb is the blast wave velocity and «hii 
and n sp are the densities of the H II region and the protru- 
sion, respectively. The fastest transmitted shocks are expected 
for head-on collisions, while slower velocities occur for tan- 
gential shocks. The shock velocities therefore depend on both 
the density and geometry of the spot. As a consequence, the 
transmitted shocks will have a wide range of velocities (~ 
10 2 - 10 3 km s~'). However, only a fraction of these shocks 
will contribute to the observed UV/optical emission from the 
spot dPun et al. 1120021) . 



Since 1997 an in creasing number o f bright spots have been 



observed in the ER (Suaermanet al. 2002). This has marked 



the transition between the SN state and the SN remnant state, 
and for the first time we have the opportunity to study such 
a transition. A central issue for understanding the physics of 
the interaction is the hydrodynamics and physical conditions 
of the shocked region. High spectral resolution in combination 
with good spatial resolution is here invaluable. To address these 
issues we present in this paper high resolution optical spec- 
troscopic VLT/UVES data, taken at 2002 October 4-7 (days 
5702-5705). In a previous paper we have discussed some spe- 
cial aspects of these data, in particul ar the presence of a numbe r 



of coronal lines from [Fe X-XIV] ( Groningsson et al. 2006). 



Here we discuss the observational analysis in greater detail, 
and in particular we concentrate on an analysis of lines from the 
low and medium ionization stages, as well as the narrow line 
component. In a forthcoming paper we will discuss the time 
evolution of the line emission using data from several epochs. 

In Sect. 2 we discuss the observations and data reductions. 
In Sect. 3 we describe how to separate out emission lines both 
from different regions such as unshocked and shocked gas and 
different spatial positions of the ER. The analysis of these data 
is done in Sect. 4 where we also describe what conclusions 
can be drawn from the line profiles and nebular analysis. The 
results are discussed and summarized in Sect. 5. 

2. Observations 

Service mode observations of supernova remnant (SNR) 
1987 A were obtained on 2002 October 4-7 using the cross- 



dispersed Ultraviolet and Visual Echelle Spectrograph (UVES) 
at the ESO Very Large TelescopeQ at Paranal, Chile. Details of 
the observations are given in Table Q] The spectrograph sepa- 
rates the light beam into two different arms. One arm covering 
the shorter wavelengths (3030 - 4990 A) has a CCD detector 
with a spatial resolution of 07246/pix. The arm covering the 
longer wavelengths (4760 - 10 600 A) has two CCDs with spa- 
tial resolutions of 0'.' 182/pix. The UVES spectrograph was used 
with two different dichroic settings and hence covering in total 
the wavelength range between 3030 - 10600 A with gaps at 
5770 - 5830 A and 8540 - 8650 A. 

To encapsulate Spot 1 and to minimize the influence from 
the bright nearby stars a 078 wide slit was centered on the SN 
and rotated to PA = 30°. For these reasons and to be able to 
separate the north and south parts of the ER we requested an 
atmospheric seeing better than 078 for the observations (see 
Fig. [TJ. Since there are no real point sources inside the slit it 
is difficult to get a direct measure of the image quality from 
the data. Instead we chose to use the DIMM as an estimate. 
It typically overestimates the actual seeing by ~ 0715 but can 
serve as an upper limit. From Table Q] we see that the average 
DIMM seeing for these observations was ~ 078. After a careful 
check we concluded that the data taken also at somewhat poorer 
DIMM seeing were good enough for our purposes, which are 
mainly to separate the north and south parts of the ER. (See 
also Appendix 1 for a discussion on the influence of seeing.) 

For the slit width used for these observations the resolving 
power is ~ 50 000, which corresponds to a spectral resolution 
of ~ 6 km s _1 . 

2.1. Data reduction 

The calibration data, according to the UVES calibration plan, 
include bias, wavelength calibration spectra, flat fields and 
spectrophotometric standard stars. Using the UVES pipeline 
version 2.0 as implemented in MIDAS, we produced the cal- 
ibration sets needed for bias subtraction, flat-fielding, order ex- 
traction and wavelength calibration of the science data. 

In order to produce reduced wavelength calibrated 2D spec- 
tra, the UVES pipeline command REDUCE/SPAT was used. 
The orders of the reduced and calibrated 2D spectra suffer 
from high noise levels near their ends, and especially so at 
the blue end. Moreover, by comparing consecutive orders of 
standard stars, we found that the flux levels at the overlapping 
regions differed significantly from each other. That could of 
course affect the relative fluxes of lines located near the ends 
of the orders. Since the automatic merging of the orders within 
REDUCE/SPAT did not give satisfactory results we chose to do 
the merging outside the pipeline. We first excluded both ends 
of each order and then applied a linear weighted averaging be- 
tween the remaining overlapping regions for every order of the 
2D frames before merging them. The overlap after cutting the 
ends of the orders were ~ 10 A except for reddest part of the 
spectrum where the overlap was too small for making adequate 
mergings. Hence, those regions were removed from the final 
ID spectrum (see Fig. [3]) The accuracy in the wavelength cali- 
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Fig. 1. Upper panel: HST/WFPC2 [N II] ,16583 image of the 
triple ring system of SN 1987A from 2002 May 10. Overlaid 
is the slit position of our VLT/UVES observations (PA = 30°) 
from 2002 October. The slit width is 0V8 and the image size 
8'.'0 x 8'.'0. Lower panel: Flux of [N II] ,16583 along the slit in 
our 2002 October spectrum. The north (left) and south (right) 
parts of the inner ring result in two prominent peaks. Note also 
the two outer ring components (marked by "OR") on either side 
of the inner ring. 



bration for the merged 2D spectra was checked against strong 
night sky emission lines and the systematic error was found to 
be ~ 0.01 A over the whole spectral range, which is also guar- 
anteed by the UVES tearrH. In order to transform the wave- 
length calibrated data to a heliocentric frame we used the task 
RVCORRECT within the IRAF package to calculate the ra- 
dial velocity component of the instrument with respect to SN 
1987 A due to the motion of the earth. 

In addition to the sky emission lines the resulting 2D 
spectra show strong background emission from the LMC. 
In order to remove this source we used the IRAF routine 
BACKGROUND. We concluded that a simple linear interpola- 
tion could not sufficiently reduce the LMC background. Instead 
we chose a second order Legendre polynomial interpolation as 
an estimate of the background level. This interpolation proved 



http://www.eso.org/instruments/uves/ 



to be a good estimate for the background emission in most 
cases. However, for [O II] and [O III] where the background 
is strong, the procedure was not able to remove all the back- 
ground emission and some background artifacts in form of 
weak emission features may still remain after the subtraction. 
Nevertheless, since the background emission component has a 
different systemic velocity than the SN most of it is spectrally 
resolved from the ring emission, and does not significantly af- 
fect our results. 

To extract ID spectra from the 2D spectra we fitted the 
double peaked spatial profile (see Fig. [T) with the sum of two 
Gaussians. For the best fit, i.e. the least-square solution, each 
Gaussian represents one part of the ER i.e. north or south. The 
two-gaussian model was then extracted in the dispersion di- 
rection. The resulting ID spectra were average-combined with 
weights proportional to their exposure times. Within this pro- 
cess we applied a sigma clipping routine which effectively re- 
jected outliers of the data such as cosmic rays which heavily 
contaminated the individual spectra. 

The flux calibration of the data was performed by using the 
UVES master response curves and the atmospheric extinction 
curve provided by ESO. Different response curves for differ- 
ent standard settings and epochs are provided. To check the 
quality of the flux calibration we reduced the spectrophotomet- 
ric standard star spectra with the same calibration sets as for 
the science spectra. The width of the slit used for the standard 
star observations was lO'.'O and slit losses should therefore be 
negligible. The extracted flux-calibrated standard-star spectra 
was compared to their tabulated physical fluxes. The fluxes ob- 
tained by using the master response curves showed to be accu- 
rate to ~ 10% in level and shape in general. The science spectra 
were corrected for this discrepancy. 

The data were observed with an atmospheric dispersion 
corrector and thus we consider the slit losses to be close 
to monochromatic. However, another source of uncertainties 
could be due to possible misalignment of the relatively narrow 
slit (0'.'8) for the science observations which could cause sub- 
stantial losses in flux. Investigations by eye of the slit images, 
provided as a part of the UVES calibration plan, revealed that 
the slit was accurately placed across the SN. Nevertheless, flux 
losses are still, of course, seeing dependent. To estimate how 
sensitive fluxes are to different seeing conditions we compared 
various line fluxes of the individual science exposures and con- 
cluded that the fluxes differ by < 10% between the different ex- 
posures. From these measurements, and the results in Appendix 
1 where we have compared our UVES fluxes to data taken with 
the Hubble Space Telescope at similar epochs, we are confident 
that the systematic error of the absolute flux should be less than 
20 - 30%. As discussed above, the relative fluxes should, how- 
ever, be accurate to ~ 10-15%. 

To create the final spectrum the combined spectra from all 
settings were merged by applying a linearly weighted averag- 
ing between the overlapping regions for the settings. 

3. Results 

From our data we were able to spatially separate the northern 
and the southern parts of the ER in the reduced 2D frames (see 
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Table 1. VLT/UVES observations of SN 1987A and its rings at days 5702-5705 past explosion. The airmasses of these obser- 
vations were ~ 1.5. 



Date 


Grating Setting 


Wavelen 


gth Range 


Slit 


Seeing 


Exposure time 






(nm) 


(arcsec) 


(arcsec) 


(s) 


2002-10-04 


CD#4 860 


uuu 


1 060 

1UUU 


12.0x0.8 


0.91 ±0.06 


3600 


2002-10-04 


CD#2 437 


j f j 


- 4Q0 


10.0x0.8 


0.91 ±0.06 


3600 


2002-10-04 


CD#2 437 


373 


-499 


10.0x0.8 


1.05 ±0.20 


2160 


2002-10-04 


CD#4 860 


660- 


1060 


12.0x0.8 


1.05 ±0.20 


2160 


2002-10-05 


CD#4 860 


660- 


1060 


12.0x0.8 


0.44 ± 0.05 


3600 


2002-10-05 


CD#2 437 


373 


-499 


10.0x0.8 


0.44 ± 0.05 


3600 


2002-10-06 


CD#3 580 


476 


-684 


12.0x0.8 


0.99 ±0.13 


3600 


2002-10-06 


CD#1 346 


303 


-388 


10.0x0.8 


0.99 ±0.13 


3600 


2002-10-06 


CD#3 580 


476 


-684 


12.0x0.8 


0.94 ±0.18 


3600 


2002-10-06 


CD#1 346 


303 


-388 


10.0x0.8 


0.94 ±0.18 


3600 


2002-10-07 


CD#1 346 


303 


-388 


10.0x0.8 


0.65 ± 0.05 


3000 


2002-10-07 


CD#3 580 


476 


-684 


12.0x0.8 


0.65 ± 0.05 


3000 




50 100 150 200 250 300 350 400 450 500 550 

Velocity (km s -1 ) 

Fig. 2. The 2D background subtracted spectrum from 2002 October 6 showing Ha. The vertical direction of the plot is along the 
slit, with the northern part in the upward direction. The total spatial extent is 12 arcsec. The velocities are with respect to the 
rest wavelength of Ha. Note how the two intersections of the slit with the ER and the outer rings are clearly visible around the 
rest velocity of the supernova, as well as the weaker emission from the outer rings. The broader (fewxlO 2 km s _1 ) emission (the 
"intermediate" component) is asymmetric, with the blue side dominating for the northern ER, and the red side for the southern 
ER. 



Fig. [2]) to extract ID spectra. Figure [3] shows the full spectrum 
from the northern part of the ring, with a number of the stronger 
lines marked. 

First we note the very broad Ha line, with an extension of 
15 000 km s" 1 . As discussed bv lSmifh et"ai](l2005l) : [Heng et al 



(120061) . this most likely originates at the reverse shock, propa- 
gating back into the ejecta. Also, at ~ 7300 A there is a broad 
component most likely from [Ca II] A Al '292, 1 '324. These will 
be discussed in detail in a future paper. 

Most of the emission lines show a narrow component 



(FWHM 



10 km s ) on top of a broader (HWZI 



300 km s _1 , henceforth called the "intermediate" component) 
(see Fig. 0). We interpret the narrow component as the emis- 
sion from the unshocked CSM and the inte rmediate compo - 
nent to come from the shocked CSM (see also lPun et al. 2002). 
To separate the emission from these two components we have 
made an interpolation of the intermediate component using 
least squares spline approximations (see Fig . |4j) . Another prob- 
lem that had to be addressed was the blending of different in- 
termediate line components. In Fig. [4] we show a typical ex- 
ample of how the deblending process was done. Here [O III] 
/15007 was separated from He I /15016 by fitting and scaling He 



I /15876 to the latter line. Assuming that these two He I lines 
have similar profiles, the He I /15016 line could then easily be 
subtracted from [O III] A5007. 

The line profile of the narrow components is mainly due 
to thermal broadening at least for the lightest elements such as 
H and He. However, since the ER is an extended object, the 
observed profile will deviate from a simple Gaussian, looking 
more like a skewed Gaussian (see Fig. [5]). By contrast, the line 
profile of the intermediate component is dominated by shock 
dynamics and there is no reason for it to be Gaussian at all. 
Figure g] shows [O III] A5001 and He I 5016 for the two com- 
ponents. In Sect. 4 we will discuss the physical cause of the 
line profiles in more detail. 

Identifications of all emission lines are given in Table [3] 
and measurements of the velocities and fluxes of the most 
significant emission lines are presented in Tables |4j-[6] The 
values listed in the tables were obtained by fitting a sum of 
Gaussians to the line profiles and the background was estimated 
by a third-order polynomial interpolation. The values given in 
brackets are the statistical lcr uncertainties of these fits. We 
note, however, that the systematic errors in the relative fluxes 
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Fig. 3. Extracted ID spectrum from the northern part of the ER. The most significant/important lines are marked. The gaps in 
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Fig. 4. The intermediate and narrow components of [O 
III] ,15006.8 and He I A5016 for the northern part of the ER. 
The dashed lines show the estimated zero intensity levels for 
the narrow and intermediate components, respectively. The red 
wing of the intermediate component of [O III] is blended with 
He I A 5015.7. He I A 5875.6 is used as a template for He I A 
5015.7 and the dotted line shows the result of the deblending 
procedure. Note also [Fe II] /15020.2 that gives the excess flux 
in the red wing of He I A 5015.7. The wavelength scale is cor- 
rected for redshift. 
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Velocity (km s _1 ) 
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320 



Fig. 5. [Fe II] Al 155 line profiles from the northern (thick line) 
and southern (thin line) parts of the ER. The velocity scale is 
heliocentric. The lines appear as skewed Gaussians due to con- 
tributions from the opposite side of the ring in both cases. See 
text for more details. 



are 10-15%, which in many cases dominate over the statistical 
errors. 

The velocity ranges of the intermediate components (Vbi ue 
and V lec |) are defined as the ranges where the flux is more than 
5% of the peak flux. The reason for using this, rather than the 
maximum extent of the line, is that the latter is very sensitive 



Table 2. The dereddened B aimer line fluxes relative to HB for 
the unshocked gas. Only the lcr statistical errors are given. 
Systematic errors are discussed in the text. 



Position 


Ho- 


H/3 


Hy 


H<5 


north 


354.2 ± 1.1 


100 


47.0 ± 0.3 


22.5 ± 0.2 


south 


322.8 ± 1.0 


100 


43.5 ±0.3 


17.9 ±0.1 


CaseB 


287 


100 


46.6 


25.6 



to the signal-to-noise ratio (S/N) of the line, which varies by 
a factor of ^ 100. With this measure we can therefore better 
compare lines of very different fluxes. The maximum extent is 
discussed in Sect. 14.2.31 

The peak velocities for the narrow components are with re- 
spect to the heliocentric motion, whereas the peak velocities 
for the intermediate components are relative to the peak ve- 
locity of the corresponding narrow component. This provides 
a natural frame for the ejecta impact on the ER. For the ex- 
tinction we adopted a redde ning of E(B - V) = 0.16 mag 
dFitzpatrick & Walbornl Il990h with E(B - V) = 0.06 mag 



from the Milky way (e.g., IStavelev-Smith et al. 1 120031) and 
E(B - V) = 0.10 mag from the L MC. The reddening law was 
taken from Cardelli et al. (1989) using Ry = 3.1. 



4. Analysis 

4.1. The narrow component 

The dereddened Balmer lines ratios (relative to HB) are listed 
in Table [2j together with the theoretical v alues given for cas e 



B at 10 4 K and a gas density of 10 4 cm 3 dOsterbrocklll989l) . 
In general, the ratios we find are close to those predicted by 
Case B theory, except for our high ju a / jiyj ratios. For both the 
unshocked and shocked (cf. Sect. I4.2.U gas, the HB flux ap- 
pears to be stronger relative to the other Balmer lines on the 
southern side of the ring than the northern, maybe by ~ 10%. 
Given the systematic uncertainty in the relative fluxes, this dif- 
ference is hardly significant. We note, however, that collisional 
excitation of Ha may be im portant for the narrow lines (see 
Lundqvist & Fransson Ill996h . 

To estimate the temperature and density of the unshocked 
ER, we have adopted a nebular analysis. For O II, S II and S III 
we have used a five-level atom, and for N II, O I and O III a six- 
level atom. Ato mic data for N II, O III, and S II are the same 
as described in iMaran et al. I d2000 | ). and for O II we include d 



data from lOsterbrockl 11989) and McLaughlin & Bel l I dl993l) . 
Data for S III are from Hayes (T986) jGalavis et al. Id 1995b and 
Tayal(1997). 

In Fig. [7] we show the allowed temperature-density range 
for both parts of the ER, usi ng our six-level atom s for [N II] 



and [O III]. As discussed by IMaran et al. I (2000), the narrow 



emission lines at these late epochs are expected to mainly come 
from gas with density significantly below 10 4 crrT 3 . Below 
this density our [O III] temperatures should be in the range 
(2.5 ± 0.3) x 10 4 K, and the [N II] temperatures should be con- 
fined to the range (1.1 + 0.1) x 10 4 K, where the errors stated 
include statistical errors as well as a +10% systematic error (see 
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Table 3. Emission lines grouped by ions. 



Ion 


Wavelengths ( A) 


H I 


3691 6b 3697 2 3703 9 3712 3721 9 3734 4 3750 2 3770 6 3797 9 3835 4 3889 3970 1 4101 7 4340 5 
4861.3, 6562.8, 8438.0, 8467.3, 8502.5, 8665.0, 8750.5, 8862.8, 9014.9, 9229.0, 9546.0, 10049.4 


He I 


3187.7, 3705.0, 3819.6, 3888.7, 4026.2, 4120.8b, 4143.8, 4387.9, 4471.5, 4713.2, 4921.9, 5015.7, 5047.7, 5875.7, 
6678.2, 6855.9b, 7065.7, 7281.4, 8444.5a, 8530.9a, 8662.2a 


Hell 


3203.1b, 3710.4a, 4541.6b, 4685.7, 5411.5, 5431.8b, 6527.1a, 8236.8a 


CI 


8727.1b 


N 1 


3466.5b, 5197.9, 5200.3 


nh 


5754.5, 6547.9, 6583.3 


01 


5577.3, 6300.3, 6363.8, 7254.2-5a, 8446.3-8a 


on 


3726.0a, 3728.8a, 7318.9, 7320.0, 7329.7,7330.7 


om 


4363.2, 4958.9, 5006.8 


Ne III 


3868.8, 3967.5 


NelV 


4714.2, 4725.7b 


Ne V 


3345.9a, 3425.9 


Mgl 


4562.7a, 4571.1 


Sill 


5041.0a, 5056.0a, 6347.1a, 6371.4a 


SII 


4068.6, 4076.4, 6716.4, 6730.8 


S III 


6312.1,9068.6, 9530.6 


cm 


9123.6a 


ArUI 


7135.8, 7751.1 


ArlV 


4711.3a, 4740.1a 


Ar V 


7005.7a 


Call 


3933.7, 7291.5, 7323.9 


Ca V 


5309.1a, 6086.4a 


Cr II 


4580.8, 4581.1, 8000.1, 8125.3, 8229.7, 8308.5, 8357.6 


Fe II 


4114.5b, 4177.2b, 4244.0, 4276.8, 4287.4, 4305.9b, 4319.6b, 4346.9b, 4352.8b, 4358.4b, 4359.3a, 4372.4b, 4416.3, 
4452.1b, 4457.9b, 4639.7b, 4774.7b, 4798.3b, 4814.5, 4874.5b, 4889.6, 4905.3, 4950.7b, 4973.4b, 5072.4b, 5111.6, 
5158.8, 5220.1, 5261.6, 5268.9b, 5273.3, 5296.8b, 5333.6, 5376.5, 5433.1b, 5477.2b, 5494.3b, 5527.6, 5721.3, 5745.7, 
6440.4b, 6872.2b, 6896.2b, 6944.9b, 7155.2, 7172.0, 7388.2, 7452.5, 7637.5, 7686.9, 8891.9, 9051.9, 9226.6, 9267.6 


Fe III 


4658.1, 4701.5, 4754.7, 4769.4b, 4881.0a, 5011.3a, 5270.4a 


Fe V 


4227.2a 


Fe VI 


4972.5a, 5277.8, 5335.2a 


Fe VII 


3586.3b, 4893.4a, 5720.7, 6087.0 


FeX 


6374.5b 


FeXI 


7891.8b 


FeXIV 


5302.9b 


Nill 


6666.8b, 7377.8, 7411.6b 



a Only a narrow component is detected. 

b Only an intermediate component is detected. 



Sect. 2.1). These res ults are in general in very good agreement 
with the estimates of M aran et al. (2000), who obtained [N II] 
temperatures in the range 9.5 x 10 3 - 1.3 x 10 4 K and [O III] 
temperatures in excess of ~ 1.7 x 10 4 K, using HST/STIS data 
from 1998 November 14, with perhaps a slightly higher [N II] 
temperature in the western ER (PA = 283°) than in the east- 



ern (PA = 103°). We find no clear evidence for a temperature 
difference between the northern and southern parts of the ER. 

To get a further handle on the temperature, we can use 
[O II] X13726, 3729 in combination with [O II] ^7319-7330, 
as well as [S II] ^^4069, 4076 in combination with [S II] 
/LI6716, 6731. The results from our five-level model atoms are 
plotted in Fig. [7] As the j : i37 29 1 'j ' xhiib and j A 6i\6l jA6Ti\ intensity 



8 Groningsson et al.: Line emission from the inner ring of SN 1987A 

Table 4. Emission lines from the unshocked gas of the ER. 



North South Extinction 

Emission Rest wavel. Relative flux" V b . . Vf whm Relative flux a V b , . Vf whm correction 





( A) 






Km s ) 


Km S ^ 






^ Km s ; 


(kms- 1 ) 




[iNe VJ 


Q/IOC OA 
j4Zj.0O 


1 "7 O _i_ 

1 1 .L ± 


1 O 


OQ 1 O" _i_ 1 O 

Zol. / ± l.Z 


Zj.4j ± U. jo 


1 1 O _i_ 

lj.0 ± 


1 A 

1 .0 


on/i /i _i_ i n 
ZV4.4 ± l.U 


19.55 


± 0.57 


o nA 
Z.UO 


[ TT"1 

L<J nj 


J /ZO.Uj 


i no o _i_ 
1UV.Z ± 


1 1 

4. 1 


too n , n -i 
Zoj.U ± U.4 


1 1 01 _i_ f\ 1 

1 j.oZ ± U. 1 j 


1 m i i 
1 U 1 . 1 ± 


i n 
4.U 


on^ 3 _i_ n ^ 
ZVj.j ± U.j 


14.77 


± 0.23 


i no 

i.yo 


[ TT"1 

L<J iij 


J /Zo.oZ 


AA /I _i_ 

00.4 ± 


J. 1 


Ton 1 i A C 

ZoU. 1 ± U.j 


1/1 i n _i_ n on 
14. ly ± U.ZU 


A3 O _i_ 

Oj.o ± 


i n 

j.y 


onn a _i_ n a 
ZVU.O ± U.O 


14.87 


± 0.37 


i no 

i.yo 


[JNe ill J 


JOOO. / J 


a a n _i_ 
44. U ± 


1 

1 .0 


OO 1 1 _i_ (\ c 

Zol. j ± U.j 


1 ^ to _i_ n 1 o 
1 j. /Z ± U. lo 


O/I O _i_ 

Z4. / ± 


1 
1 .5 


ono 1 _i_ n o 
zyz. 1 ± U. / 


17.64 


± 0.28 


i n/i 

i.y4 


rc TT1 

LS 11 J 


/inAO ah 
4U0O.0U 


1 o n _i_ 
1 /.U ± 


n n 
U.V 


o o i "2 _i_ n /i 
Zol. j ± U.4 


1 i f\n _i_ n 1 < 
1 l.U/ ± U. 1 j 


1 ^ O _i_ 

1 j.Z ± 


n o 
U.o 


on 1 i _i_ n a 
ZV l.j ± U.O 


13.77 


± 0.22 


i nn 

i.yu 


rc TT1 

11 J 


/inOA K 

4U/0. jj 


j. 4 ± 


n a 
U.O 


OO 1 H _i_ n. A 

Zol. / ± U.O 


11 /ia _i_ n -in 
1 1.40 ± U. jU 


CI, 

J.l ± 


n a 
U.O 


ono o _i_ n n 

zyz.z ± u.y 


13.60 


± 0.44 


i nn 

i.yu 


Ho 


A 1 n 1 T3 

4 1 U 1 . / j 


ZU.j ± 


n o 
U.o 


oq i o" _i_ i n 
Zol. / ± l.U 


OO /I ^1 _i_ n QA 

Zo.4j ± U. j0 


1 C A _i_ 

1 j.O ± 


n a 
U.O 


onn n _i_ n o 
zvu.y ± U.o 


24.72 


± 0.29 


i on 

Ley 


Li, , 


A1AC\ A£. 

4j4U.40 


/I A A _i_ 

44.0 ± 


1 *7 


OQ 1 1 _i_ 1 C\ 

Zol. 1 ± l.U 


OA _i_ n OA 

Zj.oO ± U.ZO 


QA -1 i 

jV.4 ± 


1 A 
1 .0 


on 1 n _i_ n q 
ZV1.U ± U.o 


26.52 


± 0.24 


1 O/I 

l.o4 


TTT1 

L<J lllj 


/1QA3 O 1 

4j0j.Z1 


1 n a _i_ 
1U.0 ± 


n o 
U. / 


z/v. j ± u. / 


1 /.jj ± U.J 1 


A O _i_ 

0. / ± 


n a 
U.O 


onn o _i_ i n 
ZVU.0 ± l.U 


23.97 


± 0.64 


1 O/I 

l.o4 


I I .. TT 

He 11 


40oj. / 


1U.1 ± 


1 .4 


Zoj.U ± l.Z 


1 ^A _i_ 1 T7 

ZI.jO ± l.Z/ 


O 1 _i_ 

/.I ± 


n o 
U. / 


on^ A _i_ 1 A 

zyj.o ± i.o 


30.39 


± 1.41 


i 

1. /j 


HB 


A A 1 3 O 

4o01. JZ 


i nn n 
1UU.U 




oon n _i_ (\ c 
ZoU.V ± U.J 


oa n/i _i_ n 1 1 
Z0.U4 ± U. 1 1 


no n _i_ 
y f.y ± 


A 

Z.O 


onn a _i_ n a 
zyu.o ± U.O 


27.36 


± 0.16 


1 O 1 

1. /I 


[ f \ TTT"1 

LO lllj 


/inco n 1 
4Vjo.V 1 


O 1 3 _i_ 


o i 
Z. 1 


oo i n _i_ n a 
Zol.U ± U.4 


1 ^ An _i_ n 1 1 
1 j.oy ± U. 1 1 


jU. 1 ± 


1 7 
1 . J 


on 1 o _i_ n z 
Ly l.o ± U.j 


22.38 


± 0.14 


1 AO 

l.Oo 


TTT1 

L<J lllj 


cnnA o/i 


o/io o _i_ 
Z4/.Z ± 


j. 1 


Zol. 1 ± U.j 


1 ^ oo _i_ n nA 
1 j.oo ± U.U0 


i za n _i_ 
1j4.V ± 


3 C 
J.J 


ono c _i_ n ^ 
ZVZ.j ± U.j 


22.85 


± 0.10 


1 AO 
1.0/ 


L<J lj 


jj / /.j4 


U.j ± 


n o 
U.Z 


too A _i_ \ A 

Zoj.4 + 1.4 


1 A AA _i_ i n 

14.00 ±1.1/ 


U.J ± 


n o 
U.Z 


oon <j _i_ 1 /i 
Zoy.J ± 1.4 


15.19 


± 1.23 


1 CO 

l.J / 


r\T TT"1 

LJN 11J 


/ j4. jy 


on a _i_ 
ZU.O ± 


n a 
U.O 


oo i ^ _i_ n "3 
Zol. J ± U.J 


1 1 An _i_ n nn 

i j.oy ± u.uy 


1 n i _i_ 
IV. J ± 


n c 
U.j 


ono n _i_ n a 
ZVZ.U ± U.4 


16.21 


± 0.16 


1 C/1 

1. j4 


I I ,. T 

He 1 


JO / J.Oj 


1 O" Q _i_ 
1 /.J ± 


n a 
U.O 


oon o _i_ n. z 
ZoU.o ± U.j 


1 a o^ _i_ n 1 < 
10. oj ± U. 1 j 


1 a n _i_ 
10. y ± 


n a 
U.O 


onn i _i_ n 
ZVU. j ± U.j 


18.63 


± 0.16 


1 C3 

1. JJ 


LU lj 


A3nn ■m 
OjUU. jU 


/i o n _i_ 
4o.U ± 




OO 1 A _i_ r\ A 

Zol. 4 ± U.4 


1 1 ao _i_ n 1 o 
1 1.0/ ± U. IZ 


An o _i_ 
4/. / ± 


i i 
1 .4 


on 1 o _t_ n /i 
zyi.z ± U.4 


14.12 


± 0.16 


1 A O 

1.4/ 


|> lllj 


A3 1 O HA 

Oj 1Z.U0 


1 A _i_ 

z.o ± 


U.J 


oon a _i_ n n 
ZoU.O ± U. / 


1 o oo _i_ n a o 
IZ. /o ± U.4/ 


1 o , 

l.J ± 


n o 
U.Z 


on 1 n _i_ i i 
Ly l.y ± l.J 


18.49 


± 0.85 


1 A O 

1.4/ 


LU lj 


AQA3 

OjOj. /o 


1 A 1 _i_ 

10.1 ± 


n a 
U.O 


OO 1 1 _i_ f\ 1 

Zol. 1 ± U.J 


11 Af\ _i_ n nn 
1 1.4U ± u.uy 


1 A O _i_ 

10. / ± 


n i 
U.4 


on 1 a _i_ n /i 
zyi.o ± U.4 


14.50 


± 0.12 


1 A O 

1.4/ 


r\T TT"1 

LJN 11J 


0j4o.Uj 


3^n n _i_ 
jjV.V ± 


"7 1 


oq i n _i_ n i 
Zol.V ± U.J 


1 o _i_ n fi? 

iz.yj ± u.uj 


j jZ. / ± 


A o 

o.y 


ono c _i_ n i 
zyz.j ± U.j 


14.88 


± 0.05 


1 A C 

1.4j 


tia 


acao on 

OjOZ.oU 


A 1 C\ 1 _i_ 

4iy.z ± 


a n 
O.V 


toa n _i_ n a 

zou.y ± u.4 


oo An _i_ n no 
Zo.Oy ± U.U/ 


i An o _i_ 
jOV.o ± 


_ _ 

j. / 


onn z _i_ n q 
zyu.j ± U.j 


27.82 


± 0.05 


1 A C 

1.4j 


LJN 11J 


Jo j.4j 


1 1 iz. / ± 


1 "7 

lo. / 


ZoU. / ± U.Z 


1 o O/i j- n n*5 
IZ. /4 ± U.UJ 


1 HQA A -L. 
1UV0.4 ± 


1 
lO.Z 


oo Mj.no 

zy i .z ± u.z 


14.99 


± 0.04 


1 A C 

1.4j 


T T , T 

He 1 


AAOO 1 C 
00/0. 1 J 


j.j ± 


n a 
U.O 


oo i i _i_ n n 
Zol. 1 ± U. / 


1 a no _i_ n 
10. U/ ± U. jj 


A 1 _i_ 

4. / ± 


n 2 
U.j 


on 1 a _i_ n o 

zy i.o ± u. / 


18.12 


± 0.32 


1 A A 

1.44 


L^ 11 J 


AT 1 A /I /I 
/ ID. 44 


OZ.Z ± 


1 . j 


OO 1 Q _i_ n, "3 

Zol.V ± U.j 


11 A1 j- n HA 

11.01 ± u.uo 


AO A -u 

ov.o ± 


1 
1 .0 


OQO Q _i_ n "3 

zyz.y ± u.j 


13.75 


±0.07 


1.43 


[SII] 


6730.82 


103.5 ± 


2.4 


281.7 + 0.3 


11.46 ±0.06 


115.5 ± 


3.3 


292.7 ± 0.3 


13.82 


±0.08 


1.43 


[Ar III] 


7135.79 


8.3 ± 


0.4 


279.9 ± 0.4 


12.91 ±0.15 


5.4 ± 


0.4 


290.1 ±0.7 


18.11 


±0.29 


1.39 


[Fe II] 


7155.16 


4.8 ± 


0.3 


281.3 ±0.4 


9.91 ±0.13 


4.9 ± 


0.2 


293.2 ± 0.4 


13.73 


±0.28 


1.39 


[Oil] 


7319 d 


13.1 ± 


1.0 






12.3 ± 


0.7 








1.37 


[Ca II] 


7323.89 


7.2 ± 


0.3 


279.9 ± 0.4 


9.77 ±0.13 


9.8 ± 


0.3 


291.9 ±0.4 


13.39 


±0.21 


1.37 


[Oil] 


7330 e 


11.4 ± 


0.8 






9.5 ± 


0.8 








1.37 


[S III] 


9068.6 


21.5 ± 


0.7 




12.76 ±0.10 


11.2 ± 


0.4 




18.90 


±0.17 


1.24 


[S III] 


9531.10 


55.2 ± 


2.0 




12.84 ±0.12 


27.8 ± 


0.9 




17.58 


±0.22 


1.22 



a Fluxes are relative to the flux of H/3 in the Northern ER: "100" corresponds to (31.3 ± 0.5) x 10~ 16 erg s~'cnr 2 . Fluxes are not corrected 

for the extinction. The errors are only statistical, and do not include systematic uncertaities (see text). 
b The recession velocity of the peak flux. 

c E(B - V) = 0. 16 mag, with E(B - V) = 0. 10 mag for LMC and E(B - V) = 0.06 mag for the Milky Way. 
d Applies to both lines at 7318.92 A and 7319.99 A. 
e Applies to both lines at 7329.67 A and 7330.73 A. 



ratios are rather insensitive to temperature, but the other ra- 
tios, i.e., /u73 1 9-733 \ I j AAil26,3129 and j^4069,4076/7/!i6716,6731 are 

temperature sensitive, these figures can be used to determine 
temperature and density simultaneously for [O II] and [S II]. 
This has not been done previously for the ER and provides a 
mor e accurate determination as i nput for models such as those 
of lLundqvist & Fransson I (119961) . For [O II] there is an overlap 
in temperature ((1.7 + 0.5) x 10 4 K) as well as in electron den- 
sity, (2.0+0.5)xl0 3 cnr 3 , of the two parts of the ER. For [S II], 
the overall preferred range of electron density vs. temperature 
is in this case (4 + 1) x 10 3 cnr 3 and (6.5 - 1 1) x 10 3 K. The 
[S II] temperature could be slightly lower on the southern side 
than in the northern. 



Again, we may compare our results against those of 
Maran et al. I d2000l) who estimated the electron density from 
[S II]. They obtained a value of ~ 4 x 10 3 cm" 3 for both the 
western and eastern parts of the ER, i.e., very similar to ours, 
albeit for other parts of the ring. They did not estimate any 
value for the [O II] density. Our nebular analysis is completed 
by temperature estimates using [O I] and [S III] for which we 
find (10 + 2) x 10 3 K and (1.3 - 1.6) x 10 4 K, respectively. 

In summary, we estimate temperatures in the range between 
~ 6.5 x 10 3 and up to at least 2.2 x 10 4 K. Listed in order 
of increasing temperature for the ions we have analyzed, we 
find S II, O I, N II, S III, O II and O III. There is a hint of 
a somewhat lower temperature on the southern side than the 
northern, especially for [S II] and maybe also for the Balmer 
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Emission 


Rest wavel. 
(A) 


Relative flux a 


^peak 

( km s~ 


') 


( km s~ 


') 


V c 
( km s 




Extinction 
correction 11 


TNe VI 

L 11V/ 'J 


3425.86 


14.8 ± 


3.6 


-136.3 ± 


21.4 


-259.1 ± 


34.9 


202.9 ± 


43.3 


2.06 


[Ne III] 


3868.75 


111.1 ± 


3.6 


- 87.8 ± 


6.2 


-263.1 ± 


4.0 


216.9 + 


8.9 


1.94 


rs III 


4068.60 


171.5 ± 


4.6 


- 77.2 ± 


4.3 


-240.5 ± 


4.7 


180.4 ± 


5.5 


1.90 


rs in 


4076.35 


45.7 ± 


2.1 


- 76.9 ± 


6.3 


-254.9 ± 


5.3 


203.0 ± 


11.2 


1.90 




4101.73 


105.5 ± 


2.6 


- 86.8 ± 


4.7 


-261.4 ± 


2.4 


197.5 ± 


6.6 


1.89 




4340.46 


223.4 ± 


5.0 


- 88.4 ± 


4.2 


-260.0 ± 


2.5 


180.5 ± 


4.1 


1.84 


ro nn 

L w 1 1 1 1 


4363.21 


32.0 ± 


2.0 


- 67.5 ± 


7.0 


-251.5 ± 


8.0 


187.7 ± 


15.1 


1.84 


[Fe III] 


4658.05 


9.8 ± 


1.1 


-101.2 ± 


11.5 


-270.9 ± 


17.3 


204.9 ± 


37.8 


1.76 


He II 


4685.7 


37.3 ± 


1.7 


- 94.7 ± 


6.8 


-265.9 ± 


5.8 


182.8 ± 


13.1 


1.75 


H8 


4861.32 


524.4 ± 


8.2 


- 85.9 ± 


3.2 


-258.7 ± 


1.1 


189.1 ± 


1.8 


1.71 


[O III] 


4958.91 


45.7 ± 


1.7 


- 66.1 ± 


7.3 


-243.7 ± 


5.6 


179.4 ± 


6.9 


1.68 


[O III] 


5006.84 


140.0 ± 


2.7 


- 58.5 ± 


4.4 


-252.5 ± 


2.9 


212.8 ± 


3.3 


1.67 


[Fe XIV] 


5302.86 


12.2 ± 


1.0 


-132.8 ± 


11.5 


-386.7 + 


10.7 


352.3 ± 


46.7 


1.61 


[O I] 


5577.34 


16.6 ± 


1.0 


- 90.3 ± 


7.0 


-253.7 ± 


5.2 


194.4 ± 


16.9 


1.57 


[N II] 


5754.59 


223.7 ± 


3.7 


- 89.6 ± 


3.4 


-252.7 ± 


1.3 


187.3 ± 


2.0 


1.54 


He I 


5875.63 


156.1 ± 


2.6 


- 84.3 ± 


3.5 


-253.3 ± 


1.4 


191.5 ± 


2.3 


1.53 


[Fe VII] 


6087.0 


2.3 ± 


0.4 


-102.6 ± 


16.4 


-345.6 ± 


22.4 


139.9 + 


22.2 


1.50 


[O I] 


6300.30 


237.4 ± 


3.7 


- 76.9 ± 


2.0 


-229.9 ± 


0.8 


169.0 ± 


1.8 


1.47 


rs nn 

i vJ 1 1 1 1 


6312.06 


11.4 ± 


0.8 


- 92.8 ± 


9.3 


-254.1 ± 


9.1 


185.2 ± 


18.9 


1.47 


[O I] 


6363.78 


81.3 ± 


2.6 


- 76.6 ± 


4.2 


-233.0 ± 


5.2 


184.0 ± 


7.2 


1.47 


[Fe X] 


6374.51 


9.8 ± 


0.8 


-127.8 ± 


10.9 


-348.7 ± 


11.7 


156.2 ± 


20.7 


1.47 


TN III 


6548.05 


75.6 ± 


1.7 


- 80.8 ± 


4.5 


-247.4 ± 


2.7 


178.6 ± 


3.5 


1.45 


Ho- 


6562.80 


2112.2 + 


31.0 


- 82.1 ± 


2.7 


-259.8 + 


0.9 


195.4 + 


0.9 


1.45 


[N II] 


6583.45 


227.5 ± 


3.7 


- 78.2 ± 


3.2 


-244.0 ± 


1.5 


183.0 ± 


1.6 


1.45 


He I 


6678.15 


j<S.Z ± 


1 .z 


- 89.6 ± 


4.9 


-261.9 ± 


3.5 


194.6 ± 


8.0 


1.44 


[S II] 


6716.44 


4.6 ± 


0.3 


- 45.2 ± 


8.4 


-232.2 ± 


7.2 


144.8 ± 


7.1 


1.43 


[S II] 


6730.82 


9.9 ± 


0.8 


- 61.5 ± 


7.7 


-239.2 ± 


19.6 


214.0 ± 


29.0 


1.43 


[Ar V] 


7005.67 


1.2 ± 


0.3 


- 89.5 ± 


22.3 


-347.6 ± 


35.3 


183.7 ± 


34.3 


1.41 


[Ar III] 


7135.79 


18.9 ± 


0.9 


- 73.8 ± 


10.4 


-255.8 ± 


5.7 


187.9 ± 


10.6 


1.39 


[Fe II] 


7155.16 


43.3 ± 


0.9 


- 85.9 ± 


2.9 


-245.3 ± 


3.0 


170.2 ± 


4.2 


1.39 


[Fe XI] 


7891.94 


9.7 ± 


0.9 


-125.9 ± 


9.8 


-336.1 + 


27.0 


218.1 ± 


25.1 


1.32 


[S III] 


9068.6 


15.8 ± 


0.9 


- 53.4 ± 


8.3 


-241.0 ± 


12.3 


138.6 ± 


6.6 


1.24 


[S III] 


9531.10 


57.4 ± 


1.7 


- 81.6 ± 


6.9 


-267.2 ± 


4.5 


172.3 ± 


6.1 


1.22 



a Fluxes are relative to the flux of HB of the unshocked component in the Northern ER: "100" corresponds to (31.3 ±0.5) x 10~' 6 erg s~'cirr 2 . 

Fluxes are not corrected for the extinction. The errors are only statistical, and do not include systematic uncertaities (see text). 
b The velocities are relative to the unshocked, narrow component of the emission lines. 
c The velocity where fluxes are 5% of the peak flux. 

d E(B - V) = 0. 16 mag, with E(B - V) = 0. 10 mag for LMC and E(B - V) = 0.06 mag for the Milky Way. 



lines. In a similar way we estimate densities in the range be- 



tween 



1.5 x 10 3 cm" 3 and ~ 5.0 x 10 3 cm 3 based on [O II] 



and [S II], respectively. This probably also brackets the elec- 
tron densities for the other ions (Fig. [7]). The exception is the 
[O III] emission, which models show mainly to come from gas 
with an electron density which is lower than ~ 1.5 x 10 3 crrT 



(e.g. lLundqvist & Fransson 119961) . These models show that the 
gas with low density remain in higher ionization stages longer. 
There is no obvious difference in density between the two sides 
of the ring. 



Another way of deriving the temperature of the emitting 
gas is from the width of the line profiles, where the FWHM 



line width due to thermal broadening relates to the temperature 

as 



= 21.40 



\10 4 A/ 



km s 



(1) 



Here A is the atomic weight of the emitting ion. Because the ER 
is an extended source, the UVES slit samples different parts 
of it, each of which having its own macroscopic and thermal 
broadening. Equation (HJ can therefore only provide an upper 
limit to the actual temperature. Because the thermal contribu- 
tion decreases with increasing ion mass, the line profiles of 
the heaviest ions should have minimum thermal contribution. 
We have chosen two such lines, but with different level of ion- 
ization to test the importance of spatial structure (cf. above), 
namely [Ar III] /T7136 and [Fe II] /17155. These lines also have 
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Table 6. Emission lines from the shocked gas of the southern ER. 



Emission 


Rest wavel. 
(A) 


Relative flux a 


V b 
(kms- 1 ) 


yv 
(kms- 1 ) 


( km s 


- 1 ) 


Extinction 
correction 11 


[Ne VI 


3425.86 


6.0 ± 1.7 


9.8 ± 28.6 


-247.6 ± 47.2 


260.7 ± 


46.1 


2.06 


[Ne III] 


3868.75 


29.5 ± 1.5 


31.7 ± 8.2 


-173.2 ± 6.7 


236.4 ± 


6.7 


1.94 


rs in 


4068.60 


41.2 ± 2.1 


31.7 ± 8.2 


-114.5 ± 6.6 


192.6 ± 


8.1 


1.90 


rs in 


4076.35 


12.6 ± 0.9 


32.8 ± 9.4 


-168.2 ± 9.0 


237.2 ± 


9.0 


1.90 


US 


4101.73 


26.5 ± 1.9 


18.8 ± 8.0 


— 118.6 ± 9.1 


215.2 ± 


12.1 


1.89 


Hy 


4340.46 


57.4 ± 2.0 


28.7 ± 7.2 


-147.5 ± 5.5 


210.5 ± 


5.9 


1.84 


ro mi 


4363.21 


10.4 ± 1.3 


32.9 ± 8.6 


-168.9 ± 7.8 


222.7 ± 


7.8 


1.84 


[Fe III] 


4658.05 


3.1 +0.6 


12.8 ± 18.7 


-232.3 ± 28.6 


259.6 ± 


28.1 


1.76 


He II 


4685.7 


10.1 ± 1.2 


39.5 ± 15.6 


-170.9 ± 20.0 


224.2 ± 


20.0 


1.75 


H8 


4861.32 


149.8 ± 2.9 


34.1 ± 3.7 


-175.1 ± 3.5 


218.5 ± 


2.0 


1.71 


[O III] 


4958.91 


18.3 ±0.8 


14.3 ± 6.2 


-226.3 ± 12.9 


204.0 ± 


6.0 


1.68 


[O III] 


5006.84 


55.3 ± 1.5 


16.4 ± 7.0 


-214.4 ± 5.8 


215.6 ± 


5.0 


1.67 


ro n 


5577.34 


3.5 ± 0.4 


38.5 ± 9.3 


-141.2 ± 10.0 


214.8 ± 


10.0 


1.57 


[N II] 


5754.59 


63.5 ± 1.4 


40.7 ± 5.9 


-188.0 ± 3.6 


207.9 ± 


3.0 


1.54 


He I 


5875.63 


49.7 ± 1.2 


38.2 ± 5.4 


-215.0 ± 3.0 


231.6 ± 


4.9 


1.53 


ro n 


6300.30 


69.1 ± 1.3 


32.5 ± 5.1 


-162.4 ± 1.8 


193.7 ± 


3.8 


1.47 


[S III] 


6312.06 


2.6 + 0.4 


24.7 ± 12.5 


-158.5 ± 16.8 


213.0 ± 


16.8 


1.47 


[O I] 


6363.78 


21.7 ±0.8 


27.1 ± 9.2 


-143.7 ± 5.4 


175.1 ± 


3.7 


1.47 


FN III 


6548.05 


24.3 ± 0.8 


31.5 ± 6.0 


-190.8 ± 4.1 


196.6 ± 


6.6 


1.45 


Ha 


6562.80 


582.0 ± 8.7 


36.0 ± 2.5 


-190.9 ± 1.2 


223.2 ± 


0.9 


1.45 


[N II] 


6583.45 


76.1 ± 1.4 


28.8 ± 3.6 


-184.2 ± 2.9 


194.2 ± 


1.6 


1.45 


He I 


6678.15 


10.0 ±0.5 


31.4 ± 8.1 


-188.4 ± 6.3 


252.8 ± 


6.3 


1.44 


[S II] 


6716.44 


1.3 ±0.2 


20.1 ± 10.0 


-124.3 ± 13.7 


167.1 ± 


13.7 


1.43 


[sm 


6730.82 


3.1 ±0.4 


27.5 ± 10.3 


-146.0 ± 12.4 


202.4 ± 


12.4 


1.43 


[Ar III] 


7135.79 


5.9 ±0.6 


19.6 ± 14.5 


-196.2 ±20.1 


217.6 ± 


15.7 


1.39 


[Fe II] 


7155.16 


11.2 ±0.5 


14.1 ± 6.3 


-184.2 ± 9.1 


199.8 ± 


8.4 


1.39 


[S III] 


9068.6 


6.8 ±0.8 


-0.1 ± 17.5 


-272.7 ± 24.0 


240.4 ± 


42.0 


1.24 


[S III] 


9531.10 


16.5 ± 1.4 


23.4 ± 8.5 


-275.5 ±25.7 


216.7 ± 


12.5 


1.22 



a Fluxes are relative to the flux of H/? of the unshocked component in the Northern ER: "100" corresponds to (31.3 ±0.5) x 10~' 6 erg s~'cirr 2 . 

Fluxes are not corrected for the extinction. The errors are only statistical, and do not include systematic uncertaities (see text). 
b The velocities are relative to the unshocked, narrow component of the emission lines. 
c The velocity where fluxes are 5% of the peak flux. 

d E{B - V) = 0. 16 mag, with E(B - V) = 0. 10 mag for LMC and E{B - V) = 0.06 mag for the Milky Way. 



the advantage of both being in the red part of the spectrum 
where the velocity resolution of the profiles are better than in 
the blue. To estimate the thermal contribution of the other lines 
we assumed that their macroscopic distributions over the ring 
are similar to that of the template line. The thermal component 
is then the deconvolution of the total line profile and the tem- 
plate. Table|7]shows the estimated temperatures for some lines. 



Due to their low atomic mass, lines of H and He are dom- 
inated by thermal broadening and the spatial structure is less 
important. The opposite is true for lines emitted by heavy ions 
such as sulphur. This is exemplified in Fig.[8]for Ha and [O I] 
/16300, where we in both cases have used the [Fe II] line as 
template. The observed profiles and fits to these are shown in 
the lower panels. In the upper panels we show the template 
line together with the thermal component of Ho- (top) and [O I] 
/16300 (bottom) needed to get a good fit. It is obvious that the 
error in estimated temperature becomes significant when the 
thermal width is smaller than the width of the template. As a 



matter of fact, the analysis becomes pointless for ions heavier 
than neon, and we have omitted such ions from Table [7] 

Table [7] shows that the line profiles indicate hydrogen tem- 
peratures of ~ 10 4 K, with Ha giving the highest temperature, 
regardless of template ion used and position of the ring. This 
is interesting since a temperature slightly in excess of 10 4 K is 
required for collisional excitation to come into play, as hinted 
by the Balmer line ratios. The fact that the unshocked ring is 
supposed to have a high abund ance of only partially ionize d 
hydrogen at these late stages (cf. lLundqvist & Fransson 1996), 
of course also boosts collisional excitation relative to recombi- 
nation. 

For He I we obtain similar temperatures as for the Balmer 
lines, whereas He II temperatures are in the range (2.5 + 1.0) x 
10 4 K with [Ar III] as template and about (4.0 + 1.0) x 10 4 K 
with [Fe II] as template. It appears likely that [Fe II] is a better 
template for H I and He I, whereas [Ar III], with its higher ion- 
ization potential, should be better for He II. There is an indica- 
tion that the temperature is slightly higher on the northern side 
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Fig. 6. Line profiles of the narrow emission lines from the northern (blue/thick lines) and southern (red/thin lines) parts of the 
ER. The zero velocity corresponds to the rest frame of the SN (286.5 km s~'). The left y-axis indicate the flux of the emission 
lines from the northern part of the ER and the right y-axis indicate the flux from the southern part. The line profiles have been 
smoothed using a fourth-order Savitzky-Golay filter. Note the thermal broadening of the lighter elements. 



Table 7. Derived temperatures from the widths of the narrow lines using [Ar III] ,17136 and [Fe II] /17155 as templates for the 
instrumental broadening and the expansion of the ring. 



Template: 


[Arlll] ,17136 


[Arm] ,17136 


[Fe II] ,17155 


[Fe II] ,17155 


Position: 


north 


south 


north 


south 


Emission line 


T (10 3 K) 


T (10 3 K) 


T (10 3 K) 


T (10 3 


K) 


Ho- 


10.5 ± 0.5 


8.0 ± 0.5 


14.6 ±0.6 


12.4 ± 


0.6 


H/3 


8.0+ 0.5 


6.9 ± 0.6 


11.6 ±0.6 


11.2 ± 


0.7 


Hy 


8.0+ 1.0 


6.5 ± 0.5 


12.0 ± 1.0 


10.7 ± 


0.9 


H<5 


9.0+ 1.0 


4.3 ± 1.3 


13.0 ± 1.0 


8.0 ± 


1.0 


He 1 ,15875.7 


7.0 ± 1.0 




13.3 ± 1.8 


11.5 ± 


1.5 


He 1 ,16678.2 


5.5 ± 1.5 




13.8 ± 1.8 


10.8 ± 


1.8 


He II ,14685.7 


28.0 ± 8.0 


25.0 ± 10.0 


40.0 ± 5.0 


40.0 ± 


10.0 


[N II] ,15754.5 






22.0 ± 2.0 


20.0 ± 


2.0 


[N II] ,16547.9 






16.3 ± 1.8 


8.8 ± 


1.8 


[N H] ,16583.3 






14.5 ± 1.5 


9.5 ± 


1.5 


[O I] ,16300.3 






10.0 ± 1.5 


3.5 ± 


1.5 


[O I] ,16363.8 






9.0 ± 1.5 


6.0 ± 


3.0 


[O II] ,13726.0 






25.0 ±3.0 


9.0 ± 


2.0 


[O II] ,13728.8 






29.0 ± 3.0 


11.0 ± 


2.0 


[O III] ,14363.2 


27.0+ 7.0 


55.0 ± 15.0 


51.3 ±6.3 


110.0 + 


20.0 


[O III] ,14958.9 


17.0 ± 3.0 


49.0 ± 6.0 


35.0 ±3.0 


105.0 ± 


10.0 


[O III] ,15006.8 


17.0 ± 3.0 


43.0 ± 8.0 


37.0 ± 3.0 


110.0 ± 


10.0 


[Ne III] ,13868.8 


23.0 ± 3.0 










[Ne V] ,13425.9 


40.0 ± 10.0 


15.0 ± 10.0 
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Fig. 7. Electron density vs. temperature for the dereddened line ratios of emission from the unshocked gas. The left pan- 
els show the emission from the northern part of the ring and the right panels are for the southern part. The line ratios are: 
[O I] ^6300,6364/^5577, [O II] (iA7319 - 7331/^3726,3729 and A3729//13726), [O III] ^4959,5007/^4363, [N II] 
/U6548, 6583M5755, [S II] (^4069,4076/^6716,6731 and ^6716/^6731), [S III] (X19069, 9531/^63 12). Each solution curve 
is embraced by error curves of at least 15% to account for the systematic relative flux uncertainties. The exception is for 
/13729//13726 and /16716//16731 ratios for which the systematic errors should be small. Instead these are embraced by the l<x 
statistical error curves. Note how the rather temperature insensitive ratios A3129/ A3126 and /16716//16731 together with the 
temperature sensitive ratios AA7319 - 733 1//L13726, 3729 and /L14069, 4076//L16716, 6731 can constrain both the density and 
temperature of [O II] and [S II]. 



for the gas emitting H and He lines, which is consistent with 
the Balmer line ratios (see Table |2). It is tempting to suggest 
that this higher temperature is a result of higher shock activity 
on that side of the ring, which could photoelectrically excite 
the unshocked gas. 

For the metal lines the thermal broadening is much smaller 
than for H and He and the method becomes less exact. 
Nevertheless, it can pick up trends in the estimated tempera- 
ture, and it offers an independent test of line ratio temperatures. 
Such a trend is seen for the lightest metal ion N II for which 
the northern side gives a higher temperature than the southern, 
i.e., in agreement with line ratio temperatures. It is also inter- 
esting that [N II] A5755 gives a higher temperature than the 
/L16548,6583 doublet. The latter comes from a lower excita- 



tion level and is biased to come from regions of lower temper- 
ature and density. The same is seen for the [O III] counterpart 
[O III] ,14363 compared to [O III] AA4959, 5007. Whereas the 
[N II] temperatures are in fair agreement with those from line 
ratios, the [O III] temperatures are inconsistent with the line 
ratio temperature for the southern side. This is most likely due 
to the fact that the [O III] lines are much stronger on the north- 
ern side and that there could be "leakage" of emission from the 
northern side to the blue wing of the emission from the south- 
ern side. The southern side line profile thus becomes too broad 
and the estimated temperature too high. The same is true also 
for the [Ne III] lines. Table [7] clearly shows the importance of 
using the correct template line for [O III], and we have retained 
the figures in the last two columns only to highlight these ef- 
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Table 8. The dereddened Balmer line fluxes relative to H/3 
for the shocked gas. Only the lcr statistical errors are given. 
Systematic errors are discussed in the text. 
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Fig. 8. The upper left panel shows [Fe II] /T7155 north 
(blue/thin line) and a thermal (gaussian) line profile for hydro- 
gen with a temperature of 14.6 x 10 3 K (red/thick line). The 
lower left panel shows the result after convolving [Fe II] with 
the thermal component (red/thick line) vs. the observed Ha line 
profile at the corresponding position (blue/thin line). The right 
panels show the same for [O I] /16300, and a thermal line profile 
with a temperature of 10 x 10 3 K for oxygen. The zero velocity 
corresponds to the peak flux of the lines. 



fects. The only [O III] line profile temperatures to be trusted 
are for [Ar III] as template and for the northern side. Those 
temperatures are also in general agreement with line ratio tem- 
peratures. 

For the other metal lines we note that the [Ne III] /13869 
agrees with the [O III] temperatures, as expected, and that 
[Ne V] indicates a higher temperature. The uncertainty of the 
latter is, however, very large due the noisiness of that line. For 
[O I] and [O II] we see the same effect as for [N II], i.e., the 
northern side appears hotter. The suggested temperature differ- 
ence between the two sides of the ring is, however, too pro- 
nounced compared to line ratio temperatures for [O II], even 
if one would consider a difference in density of the [O II]- 
emitting gas on the two sides of the ring. It should be kept in 
mind that [O II] AJ3126, 3729 lie far out in the blue where the 
spectral resolution is worse and the line profile temperatures 
less accurate. In addition, the line ratio temperature included 
the [O II] lines between 73 1 9-732 1 A which also affects a direct 
comparison between line ratio and line profile temperatures for 
[Oil]. 

4.2. The intermediate component 

In Tables [5] and [6] we give fluxes of the intermediate line com- 
ponents. The table also provides velocities of the peak of the 
line profiles, as well as maximum red and blue velocities of 
the lines. All velocities are with respect to the corresponding 
narrow line component. 

In Figs.l9land[T0lwe show a selection of line profiles of dif- 
ferent ionization stages from the intermediate component. To 



Position 


Ho- 




w 


h 7 


HS 


north 


342.3 ± 


0.4 


100 


45.13 ±0.09 


22.71 ±0.07 


south 


328.6 ± 


0.5 


100 


41.4 ±0.2 


19.0 ±0.2 


CaseB 


287 




100 


46.6 


25.6 



facilitate a comparison with other lines we have in each fig- 
ure also included the line profile of the high S/N line of [O 
III] /15006.8. From the figures it is apparent that the S/N varies 
by a large factor from the strongest lines to the weaker, which 
have fluxes less than a percent of the strongest. We also see that 
both the profile and the peak velocities differ strongly between 
the northern and southern components (see Sect. 14.2. 3b . In gen- 
eral, the northern component has the highest flux, and we will 
therefore in the following concentrate our discussion on this 
part. 

4.2.1. Nebular diagnostics 

As for the narrow lines in Sect. I4.ll we have used the fluxes 
in Tables [5] and [6] to probe the temperatures and the electron 
densities of the shocked gas. The results from our multilevel 
model atoms are plotted in Fig. [T2] The figure panels also in- 
clude shaded areas which encapsulate likely isobars (expressed 
in pressure units of 10 10 K cm -3 ). It is, as we will discuss else- 
where, a good assumption to assume constant pressure for the 
radiative shocks propagating into the ER. The reason we find 
these isobars more likely than others is that an [O III] tem- 
perature below the ~ 10 4 K typical of photoionized gas seems 
unlikely (marking the right boundary of shaded are), and that 
the [O III] temperature should be higher than the temperatures 
of [N II] and [S II], which roughly coincides with an [O III] 
temperature of ~ 4 x 10 4 K (giving the left boundary). The 
latter temperature is more typical of shock excitation (see next 
section). From Fig.Q~2]this temperature interval for [O III] con- 
strains the electron density in the [O III] emitting region to be 
between 10 6 - 10 7 cirT 3 for both the northern and southern parts 
oftheER. 

Furthermore, Fig. [12] shows a very reasonable temperature 
structure of the emitting gas in the shaded area, with the tem- 
perature decreasing, in order, through the [O III], [N II]/[S III], 
[O I], [S II] and [N I] regions. In the [N I] region it could 
be as low as (4 - 5) x 10 3 K, although the exact temperature 
depends on the degree of ionization there. For the unshocked 
ER we found the sequence [O III], [S III], [N II], [S II] and 
[O I], which would be the same also for the shocked gas if we 
limit ourselves to the upper densities of the shaded area in Fig. 
[12] where the [S III] temperature is higher than that of [N II]. 
However, this can only be regarded as a hint since there is no 
clear reason why the temperature stratification should be ex- 
actly the same in the unshocked and shocked gas. 

In gen eral our results ag ree with the nebular analysis per- 
formed by IPun et al. I (12002), as they found roughly the same 
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Fig. 9. Line profiles of the broad emission lines (blue/thick lines) vs. [O III] A5QQ1 (red/thin lines) from the northern part of the 
ER. The zero velocities are relative to the peak velocities of the unshocked, narrow component of the emission lines. The flux of 
[O III] is indicated by the right y-axis. The line profiles have been smoothed using a fourth-order Savitzky-Golay filter. Note the 
difference in profiles between the lower ionization lines and the coronal lines. 



(within errors) elemental sequence of decreasing temperature. 
However, their data suffer from comparatively large uncertain- 
ties for some of the line ratios, making a direct comparison 
rather difficult. Nevertheless, we find good agreements for the 
[N II] and [O III] densities and temperatures. From our nebular 
analysis we find, however, con sistently higher te mperatures for 
[O I] than those estimated by iPun et al. (2002), although the 
flux ratios are similar. The origin of this is not clear. 



4.2.2. Origin of the optical lines 

A detailed discussion of the charact eristics of t he em ission 
from the shocked gas can be found in [Pun et al. I J2002I) . Here 
we only discuss the points relevant to our observations. 

It is easy to exclude thermal broadening as the cause of the 
intermediate components. A Vfwhm ~ 200 km s would in- 
dicate a temperature of T ~ 10 6 K, and at this temperature 
the observed ions would be collisionally ionized. Instead, these 
optical emission lines of the low ionization species originate 
from the photoionization z one behind the ra diative shock, with 
a temperature of ~ 10 4 K dPun et al. I l2002l) . It is therefore ob- 
vious that the shape of the line profiles must be dominated by 
the shock dynamics, rather than thermal broadening. 



We have calculated models of radiative shocks for several 
shock velocities between 100 - 500 km s _1 , using a shock 
code based on the most recent atomic data. The details of this, 
as well as mo dels for the shock em ission in SN 1987 A, will be 
discussed in Fransson et al. I d2007h . Here we only make some 
general remarks about the origin of the different lines. As an 
example, we take the V s = 300 km s _I shock model with a 
pre-shock density of 10 4 crrT 3 , similar to that of the ER. 

Immediately behind the shock the temperature for this 
shock velocit y is ~ 1.5 x 10 6 K. This i s where part of the soft 
X-ray lines (Zhek ov et al. 2005 2006) and the coronal lines 



discussed in Groningsson et al.l ( 2006) arise. For shock veloc- 
ities in the range V s = 250 - 400 km s _1 the thickness of this 
collisionally ionized, hot region is 



D, 



cool 



\10 4 cm- 3 / 1 300 km s- 1 / 



cm, 



(2) 



where «o is the electron density in the pre-shock gas, i.e., in 
the ER. At this point the temperature is ~ 10 s K. The gas then 
undergoes a thermal instability and cools to ~ 10 4 K in a very 
thin region. After this the gas temperature is stabilized by pho- 
toionization heating, and a partially ionized zone takes over. In 
this zone the degree of ionization, x e , and temperature decrease 
from x e ~ 1.0 and T e ~ 10 4 K, respectively, to very low values. 



Groningsson et al.: Line emission from the inner ring of SN 1987 A 



15 




-400-200 200 400 
Velocity (km s~ 1 ) 




-400-200 200 400 

Velocity (km s~ 1 ) 




-400-200 200 400 
Velocity (km s~ 1 ) 




0.25 



-400-200 200 400 
Velocity (km s~ 1 ) 



Fig. 10. Line profiles of the broad emission lines (blue/thick lines) vs. [O III] /15007 (red/thin lines) from the southern part of the 
ER. The zero velocities are relative to the peak velocities of the unshocked, narrow component of the emission lines. The flux of 
[O III] is indicated by the right y-axis. The line profiles have been smoothed using a fourth-order Savitzky-Golay filter. Because 
of the low flux we do not show the coronal lines for this region. 



The models also show that the time for a shock to become 
radiative, f coo i, is 

years. 



(3) 



fcool * 6 - 7 (l0 4 cm -a) (300km s-i) 

In the absence of a magnetic field the density in the pho 
toionized region is 



« P h 



10' 



UOOkms- 1 / \ 



10 4 K 



(4) 



where T p h is the temperature in the pho toionized zone. If there 
is a sufficiently strong magnetic field this may limit the com- 
pression behind the shock to ~ (87rpo) 1/2 V s /.Bo, where po and 
Bo are the d ensity and magnetic field in the pre-shock gas, 
respectively (iCox 
10 



1972). For V s 



^ 4 cm" 3 



300 km s 

we therefore expect a density of ~ 10 7 cm 



1 and «o ~ 
3 which 



is in good agreement with the density we infer in this region, 
4 x 10 6 - 2 x 10 7 cm" 3 fSect. BXfl , There is therefore not from 
this observation any need for a dominant magnetic field. 

Besides the coronal lines, also the [Ne V] lines come from 
the cooling, collisionally ionized region with (1 - 3) x 10 s K. 
In terms of emissivity, [Ne V] takes over roughly where [Fe X] 
drops. The [O III] and [Ne III] lines arise from a wide range 
of temperatures, (2 - 10) x 10 4 K, in the collisionally ionized 
zone, but also have a contribution from the photoionized zone 
at ~ 10 4 K. 

Both the H I and He I lines are dominated by the photoion- 
ized region just behind the cooling region at ~ 10 4 K. Most 
of the emission comes from the region where x e > 0.1, with a 
thickness of only ~ 10 12 cm. The low ionization lines from [O 
I], [N II], [S II], and [Fe II] all come from the photoionized 



zone, but the emission extends to a lower temperature and ion- 
ization compared to the H I and He I lines. With the exception 
of the coronal lines, the optical line emission therefore require 
radiative shocks. 

4.2.3. Line profiles 

The interpretation of the shape of the line profile is not straight- 
forward, although it is clear that it is dominated by the fluid mo- 
tion behind the shocks, while thermal broadening should only 
have a minor influence. As discussed in Sect. Q] it is, however, 
likely that there is a range of shock velocities, determined by 
the blast-wave velocity, density and the geometry of the blobs. 

Furthermore, the line profiles only reflect the projected 
shock velocities along the line of sight, and consequently, the 
FWZI velocity is only a lower limit to the maximum veloci- 
ties of the radiative shocks driven into the protrusions. At this 
epoch (2002 October) the shocks have had about ~ 5 years 
to cool since first impact, and from the FWZI of the line pro- 
file of Ha (see Tables [5] and [6]), we estimate that shocks with 
velocities < 260 km s _1 have had enough time to cool. This ve- 
locity depends on the density, V coo i oc « ' 29 . As faster shocks 
become radiative the maximum shock velocity resulting in ra- 
diative shocks is expected to increase with time. Therefore, the 
widths of the emission lines are also expected to increase. 

If we compare the line profiles of the different lines we find 
some very interesting differences. The low ionization ions such 
as Ha, He I, [O III] and up to [Ne V] and [Fe VII], all have very 
similar profiles, with peak velocities -60 to - 80 km s _1 , with 

Vbiue 260 km s _I for the northern part of the ER (Fig. 0. 

The coronal lines have, however, a considerably higher peak 
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velocity of 130 km s" 1 , and extend for [Fe XIV] to Vbi ue ~ 

-387 km s" 1 . This indicates that these two groups of lines do 
not arise from exactly the same regions. 

An advantage with this comparatively early epoch is that 
the emission from the northern part of the ER is still domi- 
nated by Spot 1. At more recent epochs, a large number of 
blobs along the ER, all with different line of sight projec- 
tions, contribute to the emission. There should therefore at the 
epoch discussed in this paper not be any serious contamina- 
tion from other spots with different line-of-sight angles and 
shock velocities. The line profile is therefore dominated by 
the geometry and density distribution around the blob. Further 
support for this comes from more rec ent observations in 2005 
with the SINFONI instrument at VLT dKiaer et al.l2007l) . These 
adaptive optics observations have lower spectral resolution, but 
higher spatial resolution than UVES, ~ 0'.'2. These observa- 
tions show that even in 2005, when many blobs contribute to 
the emission in the ER, there is an agreement between the peak 
velocity of the lines in the direction of the UVES slit, as mea- 
sured by SINFONI, and that obtained from the UVES observa- 
tions. In contrast to our UVES observations, the SINFONI ob- 
servations, do not show the full distribution of velocities, only 
the average at each point. 

In this context we note the difference between the north- 
ern and southern line profiles shown in Figs. 9 and 10 (see 
also table 4 and 5). While the peak velocities, V pea k, are 
-50 - (-90) km s _1 for the low and intermediate ionization 
lines from the northern part of the ER, they are in the range 
10 - 40 km s" 1 for the southern part. As discussed in Kjaer et 
al., these differences are consistent with the general geometry 
of the ER. Because of the low fluxes we can, however, not say 
much about the profiles of the coronal lines from the southern 
part. 

When we compare the line profiles of the low ionization 
lines (e.g. [O I], [S II] and [Fe II]) with the intermediate ioniza- 
tion lines (e.g. [O III] and [S III]) we find that for the northern 
part the blue wings are very similiar. However, the red wings 
of the low ionization lines are significantly weaker compared 
to the intermediate ionization lines. The opposite is true for the 
southern part. This indicates that the excitation conditions are 
dif ferent in the two r egions . 



Pun et al. I (I2002h find that a simple spherical geometry well 
explains the basic features of the line profile from Spot 1. 
They, however, find it difficult to fit the wing with velocities 
S= -150 km s" 1 with this model, and therefore invoke a 'turbu- 
lent' smoothing, possibly originating from instabilities in the 
shock structure. While this may be reasonable, it is difficult 
to understand how this can produce velocities higher than the 
shock velocity. Instead, we believe that it is more reasonable 
that the wings come from the highest velocity shocks, where 
only a small fraction of the gas, presumably that with the high- 
est densities, has had time to cool. There should then also be 
a contribution from adiabatic sh ocks, with coolin g time longer 



than that indicated from Eq. Pun et al. I d2002l) suggest that 



these could be responsible for some of the X-ray emission ob- 
served. 

This picture is supported by our observations of the larger 
extent of the coronal lines, up to ~ 400 km s _1 , compared to 




600 



Velocity (km s 



Fig. 11. The narrow and intermediate components of Ha for the 
northern part of the ER. The dashed lines show the estimated 
zero intensity levels for the narrow and intermediate compo- 
nents respectively. The blue wing of the intermediate compo- 
nent is only barely separated from the emission of [N II] /16548. 
The zero velocity corresponds to the peak flux of the narrow 
component. 



the lower ionization lines. Some of the coronal emission may 
therefore originate from adiabatic shocks, which are also likely 
to give rise to some of the so ft X-ray emission seen by Chandra 



(Zhe kov et al. 2005, 2006). If this is correct one expects the 



maximum velocity of both the coronal lines, and especially the 
low ionization lines, to increase with time. 

We note that the shock velocity should be ~ 4/3 x the 
gas velocity immediately behind the shock, as reflected in the 
line widths of the coronal lines. The velocity of the cool gas 
seen in the low ionization lines should, on the other hand, 
be close to that of the shock velocity. Therefore an extent of 
> 400 km s" 1 of the coronal lines indicates a shock velocity 
of at least > 510 km s" 1 . 

Although the emission from the coronal lines clearly comes 
from higher velocities than the lower ionization stages, it is 
important to realize that there are also in these evidence for 
gas at velocities considerably higher than Vbiue or V reQ s which 
only gave the extension to 5% of the peak flux. In particular, 
the Ha line, with the best S/N of all lines, shows clear ev- 
idence of an extension to considerably higher velocity than 
Vbiue ~ -260 km s . To illustrate this we show in Fig. QT| 
Ha on a logarithmic flux scale, From this we see that this line 

has a weak blue extension up to 450 km s" 1 , while the red 

extension is up to ~ 350 km s _1 . Therefore, at least a small 
fraction of these high velocity shocks have time to cool. Most 
likely this comes from shocked gas with densities considerably 
higher than that typical in the ER, ~ 10 4 cm" 3 . It is in this con- 
nection interesting that lLundqvist & Fransson i 19961) find evi- 
dence for densities > 3 x 10 4 cm" 3 in the ER. From Eq. (0 we 
find that for a shock velocity of 450 km s" 1 gas with a density 
of > 4 x 10 4 cm" 3 should have time to cool. 
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We finally remark that it is well known that radiative shocks 
in this velocity range are unstable to i nstabilities with wave- 
length of the order of the cool i ng len g th Jchevalier & Imamural 



19821 Strickland & Blondin 1995; SutherlandetalJ 2003) 



These may add further to the already complex line profiles. 



5. Summary 

Our UVES observations show the dynamics of the shock inter- 
action in SN 1987 A with unprecedented spectral resolution. We 
find evidence for shock velocities up to ~ 500 km s _1 . From the 
larger extent of the line profiles of the coronal [Fe X-XIV] lines 
compared to the low ionization lines, we argue that the high- 
est velocities come from adiabatic shocks, which have not had 
time yet to cool. Shocks with velocity less than ~ 250 km s _1 
are radiative and give rise to the rest of the optical lines. We do, 
however, also find evidence for a minor component of radia- 
tive shocks with higher velocity. As more and more gas cools 
we expect the width of especially the low ionization lines to 
increase. 

While the coronal lines and high ionization lines like [O III] 
and [Ne III-V] arise in the collisionally ionized hot gas behind 
the shocks, the low ionization lines come from gas photoion- 
ized by the shock radiation. The densities we derive from the 
line ratios are consistent with the large compression expected 
in a radiative shock. 

The emission at this epoch is dominated by Spot 1. In a 
subsequent paper we will discuss the evolution of these lines 
with time as more spots emerge and the shocked gas at this 
epoch has had more time to cool. 
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Appendix A: Test of flux calibration 

To test our flux calibration of the UVES spectra, we 
downloaded archival HST/STIS data for two epochs: 2000 
November 3 and 2002 October 29. The epochs are close to 
our UVES data from 2000 December 9-14 (discussed in 
Groningsson et al. 2007 in prep.) and 2002 October 4-7 (dis- 
cussed in this paper). The STIS observations are long-slit spec- 
troscopy observations made with the first-order gratings G430L 
and G750L which cover the spectral intervals 2900 - 5700 A 
(resolution 2.73 A pixeL 1 ) and 5250-10300 A (resolution 
4.92 A pixeL 1 ), respectively. The aperture was 52" x 2", and 
the slit was placed so that it covers half of the ring, with the 
dispersion direction in the north-south direction. A second set 
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Fig. 12. Electron density vs. temperature for the dereddened line ratios of emission from the shocked gas. The left panel 
shows the emission from the northern part of the ER and the right panel is for the southern part. The line ratios are: 
[O I] ^6300,6364/^5577, [O III] ^4959,5007/^4363, [N I] AA5 198,5200/^3467, [N II] ^6548,6583/^5755, [S II] 
^4069,4076/^6716,6731, [S III] ^9069,9531/^6312. Each solution curve is embraced by error curves of at least 15% 
to account for the systematic relative flux uncertainties. The shaded region marks the range for possible isobaric solutions. Four 
isobars are plotted (dashed lines) and the pressures are (in units of 10 10 K cm -3 ): 3.9(3.1), 5.0(4.1), 6.3(5.5) and 7.9(7.4). The 
values in brackets are for the southern part. 



Table A.l. Comparison of fluxes from HST/STIS and UVES data sets. 



Lines 


Grism 


Flux (ergs 


s 1 cm 2 ) 




Ratio 


Date of HST/STIS observations 






STIS 


UVES 




STIS/UVES 




Hy ,14340 


G430L 


1.06 x 10 -14 


5.74 x 10" 


ts 


1.85 


2000-11-03 


[O I] ,16300 


G750L 


1.20 x 10" 14 


5.80 x 10- 


15 


2.05 


2000-11-03 


[S H] ,1,14068,4076 


G430L 


2.35 x 10" 14 


9.83 x 10" 


15 


2.39 


2002-10-29 


H/3 ,14861 


G430L 


7.19 x 10"' 4 


2.73 x 10" 


14 


2.63 


2002-10-29 


[N II] ,15755 


G750L 


2.44 x 10" 14 


1.01 x 10- 


14 


2.41 


2002-10-29 


He 1 ,17065 


G750L 


7.59 x 10" 15 


3.56 x lo- 


15 


2.13 


2002-10-29 



of observations were made for the other half of the ring. One 
therefore obtains separate spectra for the northern and southern 
parts of the ring of SN 1987 A. These have to be added together 
to get the total ring flux. There is a marginal spatial overlap 
between the two sets of observations, but this is small enough 
(i.e., a few per cent) to be unimportant for our test. 

The pipeline-reduced STIS spectra have several cosmic 
rays, but there are fortunately "clean" lines in most parts 
of the spectrum so that we can test the flux calibration for 
the various spectral settings of UVES. The lines we chose 
in the red part of the spectrum (G750L) are [N II] ,15755, 
[O I] ,16300 and He I ,17065, and in the blue (G430L) we se- 
lected [S II] ^4068,4076, Hy and HJ3. 

To obtain the flux of a particular line in the STIS spectrum, 
we first manually cleaned the region around the line from cos- 
mic rays. We then integrated the spectrum in the dispersion di- 
rection for each spatial row that cuts across the ring. With the 
pixel scale of STIS, this means that we made spectral scans for 
up to 37 spatial rows (the exact number depends on the bright- 
ness of the ring for the line analyzed) to cover each of the full 
half ring. We then summed up the flux from all the rows in the 



spectrum and added the flux from the opposite side of the ring 
to get the total ring flux for each line. 

As the spectral resolution of STIS does not allow us to 
separate the narrow component from the intermediate-velocity 
component, the flux estimate from STIS is a sum of both these 
components. The STIS data do, however, allow us to separate 
out any underlying continuum. For the UVES spectra we there- 
fore integrated the emission in the narrow plus intermediate 
velocity component to make a direct comparison with the STIS 
flux. 

The STIS fluxes for the whole ring and ratios of STIS to 
UVES fluxes are summarized in the Table IA.1I The fact that 
the STIS flux is higher than the UVES flux is reasonable since 
the UVES slit only encapsulates part of the ring whereas the 
STIS flux is for the full ring. The ratio of STIS/UVES flux ap- 
pears to be somewhat lower for 2000 November (~ 2.0) than 
for 2002 October (~ 2.4). This could be due to temporal flux 
variations round the ring (e.g., new spots appearing between 
2000 and 2002), or may be due to a systematic error in the 
fluxing of the UVES spectra. The former is not unlikely since 
Spot 1, encapsulated by the UVES aperture, dominates the ring 
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interaction emission in the 2000 November data, whereas two 
years later, it is less dominant and various parts (also those out- 
side the UVES slit) contribute to the ring interaction emission. 
Because the STIS/UVES flux ratio does not show any obvious 
wavelength-dependent trends (see Table lA. il l, a constant multi- 
plicative factor can be used for the UVES spectrum to estimate 
the total flux from the full ring. As we have discussed, the fac- 
tor may be somewhat different for different epochs. As we here 
concentrate on our UVES data from 2002 October, we have 
chosen the factor 2.4 for this epoch. 

To check the significance of the multiplicative factor, we 
downloaded HST/WFPC2 archival data from 2002 May 10 ob- 
tained using the F656N filter. The pivot wavelength for this fil- 
ter is 6563.8 A and the band width is 53.8 A. We convolved 
the image with a 2D Gaussian function to simulate the see- 
ing of our UVES data. We tested three FWHM values of the 
Gaussian: 076, 078 and 170. The Gaussian was truncated at 5" 
as the inclusion of emission outside this radius did not change 
the results. Fi gure [ATI shows the original and a blurred image, 
simulating 078 seeing, with the 078 UVES slit overlaid on top 
of them. The total flux from the ring was assumed to be the flux 
within the larger box in the original image. (The contributions 
from the ejecta and outer rings are low in comparison with that 
of the inner ring.) We then formed a ratio of emission through 
the 078 slit and the total emission from the ring. The results are 
shown in the Table I A. 21 

For the original image (i.e., no seeing), no emission from 
the north-west and south-east parts of the ring enters the 078 
slit. With increased "seeing" emission from these parts ven- 
tures into the slit, whereas emission within the slit in the orig- 
inal image falls outside the slit. According to Table IA.2I the 
"loss" of flux measured through the slit is slightly larger than 
the "gain" since the ratio (total flux to UVES slit flux) increases 
with the degree of blurring. The effect is, however, rather small, 
with the flux within the UVES slit decreasing from 43.7% to 
40.8% with seeing increasing from 076 to 170. This means that 
seeing effects have no major impact on the absolute fluxing 
of the UVES spectra, and that spectra obtained under differ- 
ent seeing conditions can be coadded without introducing sys- 
tematic error in the fluxing. Also, since the ratios in Table lA. 2l 
show very good agreement with the ratio of 2.4 obtained from 
the comparison between UVES and STIS fluxes, we are confi- 
dent in the magnitude of this factor. The main effect of seeing 
is a redistribution of flux into and out from the area covered by 
the slit. For a time series of spectra it therefore makes sense to 
compare spectra with roughly the same average seeing to study 
changes of physical properties of the ring. 

Table A.2. Comparison of fluxes from HST/STIS and UVES 
data sets using HST/WFPC2 with F656N filter centered at the 
rest wavelength of Ha. The apertures shown on the Figure lATI 



Modeled seeing 


Ratio of aperture "2" to "1" 


unblurred 


2.31 


ff.' 6 


2.29 


a' 8 


2.35 


170 


2.45 
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Fig. A.l. Ha image of SN 1987A obtained with HST/WFPC2 in 2002 May 10 using the narrow F656N filter centered at the rest 
wavelength of Ha. The rectangular boxes are apertures which were chosen to model UVES data (marked "1") and STIS data 
(marked "2")- The image to the left is unconvolved. The other image is smoothed using 2D Gaussian functions and represent a 
seeing of 078. The lower flux cut in both images are roughly at the flux level of the outer rings in the image to the right. As the 
inner ring is nearly two orders of magnitude brighter (cf. Fig. 1) and the flux from Stars 2 and 3 fall below this flux cut at the 
position of the UVES slit, these stars have no important influence on the flux calibration. Note the dominance of Spot 1 in the 
north-east part of the ring. 



